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Abstract 



Pragmentation functions for D mesons, based on the convolution of a perturbative part, 
^ ■ related to the heavy quark perturbative showering, and a non-perturbative model for its 

hadronization into the meson, are used to describe D* production in e~^e~ and ep collisions. 
The non-perturbative part is determined by fitting the e~^e~ data taken by ARGUS and 
OPAL at 10.6 and 91.2 GeV respectively. When fitting with a non perturbative Peterson 
fragmentation function and using next-to-leading evolution for the perturbative part, we 
find an e parameter sensibly different from the one commonly used, which is instead found 
with a leading order fit. The use of this new value is shown to increase considerably the 
cross section for D* production at HERA, suggesting a possible reconciliation between 
the next-to-leading order theoretical predictions and the experimental data. 
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1 Introduction 

The study of fragmentation functions (FF) for heavy quarks has recently attracted an increased 
interest due to the large amount of data accumulated at LEP and HERA. From the theoretical 
side predictions have been obtained by combining perturbative QCD - which allows to resum 
large logarithms with a resulting milder renormalization/factorization scale sensitivity - with 
a non perturbative component which describes the hadronization of the heavy quark into the 
meson, after the perturbative cascade. 

In e~^e~ annihilation an analysis along these lines was performed by Colangelo and Nason 
'^ up to LEP energies, for both charm and beauty mesons. Due to the presence of the c(6) 
component only, their results were not applicable to the production of heavy mesons in hadronic 
collisions, where the gluon-gluon and quark-gluon scattering play an important role. Then in 
a previous analysis [Q a set of NLO fragmentation functions for D, D* mesons was given, 
including the gluon term as well, and predictions for large transverse momentum production 
cross sections were also provided. 

The aim of the present analysis is to reconsider the situation of charmed meson fragmen- 
tation functions both in e~^e~ annihilation and in photoproduction, where new data have been 
obtained at HERA. 

On the perturbative side, we consider the full set of perturbative fragmentation functions 
(PFF's) and their mixing in the evolution. This is important as the OPAL data do indeed 
show a rise in the small x region, due to the gluon splitting, which is absent in the ARGUS 
data. In addition, by parametrizing the non-perturbative component by different forms and 
fitting e~^e~ data, we study the variation of the non perturbative parameters, in particular for 
the Peterson form ||^, as related to the accompanying approximation, leading (LO) or next- 
to- leading order (NLO), used in the perturbative component. We find indeed that a NLO 
evolution favours a much smaller value of the e parameter in the non-perturbative Peterson FF 
than given in the literature. In turn this also helps reconciling the recent HERA data with the 
theoretical predictions. When however a LO evolution only is considered, as in many of the 
parton shower Monte Carlo codes used in the experimental analyses, the "conventional" value 
for e is recovered. This result can be understood by noting that the effect of parton showering, 
which is larger in a NLO analysis, softens the distribution of the partons, acting qualitatively 
as a non perturbative FF, which can henceforth behave more softly. Therefore the value of e 
used in the phenomenological analyses must be closely related to the level of the approximation 
followed in the perturbative QCD evolution. 

This paper is structured as follows: in Section |^ we recall the theoretical framework, partly 
already introduced in 0, on which this work is based. Section ^ presents the results of fits 
to ARGUS and OPAL data in e~^e~ collisions. Section ^ makes use of the non perturbative 
parameters previously determined to give predictions for D* photoproduction in ep collisions 
at HERA. Our conclusions are then given in Section ^. 



2 Theoretical Framework 

We have already introduced in Ref. [0 the theoretical framework for evaluating D mesons 
cross section within a fragmentation approach. In that paper, the following ansatz for the 
fragmentation function (FF) of a parton i into a meson D was made: 

Df(x,/i) = A'(x,/i)®D^^(x). (1) 

In this equation, D^{x,^) is the perturbative fragmentation function (PFF) for a massless 
parton to fragment, via a perturbative QCD cascade, into the massive charm quark c. D^ (x) 
is instead a non-perturbative fragmentation function, describing the transition from the heavy 
quark to the meson. Finally, the symbol ® indicates convolution, i.e. 

f(x)0g{x)^ f'-f{z)g{x/z). (2) 

Jx Z 

The formalism of PFF's has been introduced a few years ago [^], and will not be given here 
in detail. We just recall that it allows to extract from perturbative QCD (pQCD) the initial 
state conditions for the PFF's at a scale /iq of the order of the heavy quark mass m (and we 
will take Hq = m): 

D,^(x,^o)=^(l-x) + ^^^[\±:^(log^-21og(l-x)-l)] (3) 

Dl{x^^^.) = ^iif^(x^ + (1 -x)^) log 4 (4) 

Dg,g,-c{x, /io) = (5) 

where c represents here the heavy quark and g and q the gluon and light quarks respectively. 
Moreover, Cp = 4/3 and Tp = 1/2. 

The PFF's, evolved up to any scale /i via the Dokshitzer-Gribov-Lipatov-Altarelli-Parisi 
(DGLAP) equations, can be used to evaluate heavy quark cross sections in the large transverse 
momentum (p^) region (i.e. Pt ^ ttl) by convoluting them with cross section kernels for 
massless partons |§|, §, 0) subtracted in the modified minimal subtraction (MS) scheme, where 
the heavy quark is also treated as a massless active flavour and therefore also appears in the 
parton distribution functions of the colliding hadrons. This has been done in Ref. |^ for pp, 
in Ref. |^ for 7j9 and finally in Ref. [|1^] for 77 collisions. In all cases it has been shown 



how the results agree with the full massive ones (Refs. []TT|, |T^ and l]!^ respectively) in 
an intermediate px region (say from twice to four times the mass of the heavy quark). For 
larger pt they are more reliable (and hence have a smaller scale dependence) because the large 
logarithms originating from gluon emission and gluon splitting are resummed by the evolution 
of the PFF's (see Ref. [§ for a more complete discussion on this point). 

The fragmentation functions of eq. ([1|) will be also evolved with the DGLAP equations. It 
is to be noted that in doing so we assume the evolution to be entirely perturbative in character: 
we evolve the full FF's (|I|) as we would the PFF's only. The non-perturbative part of the overall 
FF's is kept fixed and determined at a given experiment. 



Indeed, the non-perturbative part of the FF's cannot be predicted by pQCD. In fact, the 
process through which a heavy quark binds to a light one to form the meson involves ex- 
changes of gluons with momenta of order Kqcd or smaller, and is therefore intrinsically non- 
perturbative. However, a few features of this function can be determined. In contrast to light 
quark hadronization, this FF is hard [14] because the meson retains a larger fraction of the 
heavy quark initial momentum. Moreover, one expects the non-perturbative FF to be squeezed 
towards x = 1 linearly in the mass of the heavy quark. This statement is proved in [^ under 
the hypothesis of softness of the hadronization process and in the infinite mass limit (see also 
T6| for a discussion on this point). 



In the following we will employ two different functional forms for the non-perturbative part 
of the fragmentation function. 

The first one is dictated mainly by its semplicity, and is given by 



Dnp{x-a,(3) = A{l-xrx'' (6) 



with 

1 /■! 



A JO 



l-x)''x'^dx = B{f3+l,a + l), (7) 



B{x,y) being the Euler Beta function. This functional form had already been employed in [|T| 
for fits to e^e~ data and was also used in our previous paper on charmed meson FF's 0. It 
is flexible enough to describe the data and has the advantage of an easily calculable Mellin 
transform, given by 



Dnp{N;a,(3) = / dxx'^ ^Dnp{x;a, (3) 
^ u 



5(/? + Ar,« + l) r(/5 + Ar)r(a + /? + 2) 



5(/? + l,a + l) r(/3 + l)r(a + /3 + iV + l)' 

with r(a;) being the Euler Gamma function. 

However, this functional form has no immediate physical motivation. A successful descrip- 
tion of e^e~ data could be not enough to ensure the correctness of the predicted cross sections 
in, say, ep production evaluated with the same non-perturbative FF, since higher moments 
could play an important role. Indeed, in e^e~ collisions it is the mean scaled energy, i.e. 
J dzzD{z) - or the second moment when talking Mellin transforms language - the most im- 
portant observable. Different FF's could therefore agree on this second moment but then have 
different higher moments which could lead to different prediction in other kinds of reactions. 

We have therefore chosen to employ also a different non-perturbative fragmentation, based 
on a physical model: the so called Peterson form [Q. It is derived by considering the transition 
amplitude for a fast moving heavy quark Q to fragment into {Qq) + q, q being a light quark. 
It reads 

Dnn{x:e) = 9, (9) 

^^ ^ x[l-l/x-e/{l-x)f ^ ' 

with the normalization factor A now given by 

1 (e2-6e + 4) f e 2-e 1 1, 1 

— = , { arctan : + arctan , : > H — m e -|- . ( 10 ) 

A (4-e)v/47^^\ ViT^l^ ViT^l^J 2 4-e ^ ' 



From the derivation one finds that the e parameter is related to the heavy quark mass by 
e ~ A^/?7i^, where A stands for a hadronic scale. Since the average scaled energy goes like 
(x) = 1 — i/e, we see it respects the prediction of scaling linearly with the heavy quark mass. 

While this form of non-perturbative fragmentation function is certainly more physical and 
the order of magnitude of its unknown parameter can be estimated from first principles, it has 
however the drawback of a much more complicated Mellin transform. The full expression is 



given in the Appendix of Ref. |T^, and will not be repeated here. 



3 Production in e+e collisions 

According to QCD factorization theorems, the cross section for the production of a hadron H 
in the e~^e~ process 

e+e-^7,Z^ifX, (11) 

at a center-of-mass energy Q = ^/s, can be written as 
da" ^^ /-i dz 



dx 



r dz { X \ 

J2 / —Ci{z, a,(/i), Q, fi)D^ [-,I^)=J2 ^^(-2' "^(/^)' Q^ Ai) ® D"i^^ A^)' (12) 

■ J X Z \ Z / ■ 



X being the energy fraction of the produced hadron, x = 2E/Q. The functions Ci{z, as{fJ.), Q, jj) 
are the so called coefficient functions, which describe the hard part of the scattering process 
and can be calculated in perturbation theory as series expansions in the strong coupling ag^jj). 
Explicit expressions up to NLO for all the coefficient functions we need can be found, for 
instance, in Ref. |l^. Since we take the partons in the hard scattering to be massless, collinear 



singularities appear, and these are subtracted in the MS scheme and reabsorbed into the 
fragmentation functions, /i is the factorization scale at which this subtraction is performed, 
which in this case we have for simplicity taken equal to the renormalization scale. The sum is 
to run on all the partons which can be considered massless in the coefficient functions. Since in 
general mass terms of the form of powers of m/Q will appear, we see that already at Q = 10 GeV 
the charm can to a good approximation be taken as massless. The same will be true also for the 
bottom quark at Q = 91 GeV, whereas its production should instead be strongly suppressed 
at the lower energy. We will therefore include four and five active flavours respectively at these 
two center-of-mass energies. 

When dealing with light hadrons the fragmentation functions can only be determined by 
comparison with experiment. Since in our case the hadron in question is instead the heavy 
meson D*, we can make use of our ansatz of eq. (|1|), and flt to the experimental data only the 
non-perturbative part of the FF's. 

We start by trying to fit the non-perturbative FF to experimental data for D*"^ production 



taken by ARGUS ||] and OPAL at 10.6 GeV and 91.2 GeV respectively. The cross section 
is evaluated by means of the formula in eq. (|T^) , the fragmentation functions are given by the 
initial conditions reported in the previous section, evolved up to the desired scale with the 
DGLAP equations to next-to-leading (NLO) order and convoluted with the non-perturbative 
component. 
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xVd.o.f 




With Sudakov resummation 




ARGUS, 


Ref. 11 0.4 4.6 




ARGUS 


0.51 ± 0.37 4.9 ± 1.7 


0.70 


OPAL 


0.30 ± 0.21 4.5 ± 1.5 


1.26 




Without Sudakov resummation 




ARGUS 


1.0 ±0.6 6.7 ±2.3 


0.86 


OPAL 


0.9 ± 0.3 6.4 ± 1.9 


1.32 



Table 1: Results for the fitting of a and /? in (1 -x)°x^ to ARGUS and OPAL 
data. Evolution is performed to NLO and with A5 = 200 MeV and /Uq 



m. 



3.1 Fits with (1 



x)"x^ 



We first perform fits with the "simple" form {l — x)°'x'^. Similar fits had already been performed 
a few years ago in Ref. [jl|. In that paper only the non-singlet component of the FF's had been 
taken into account, a valid approximation at the low energy of 10.6 GeV. When going to higher 
energy, on the other hand, the mixing with the gluons through the evolution will become more 
and more important. The OPAL data do indeed show a rise in the small x region, due to gluon 
splitting and absent in the ARGUS data. We have therefore included the full set of FF's and 
mixings in the evolution. 

As a first step, we have refitted the same ARGUS data already considered in Ref. 0]. We 
have taken A5 = 200 Mev and included in the PFF's the resummation of Sudakov terms in the 
x ~ 1 region, as described in ||^ and consistently with [|l|. A normalization factor is always 
fitted along with the parameters determining the shape of the non-perturbative FF. The results 
are shown in the upper part of Table |I|. 

They can be seen to be consistent with those obtained in Ref. [Q. It is also worth mentioning 
that the last point in the ARGUS data has not been included in our fit. In that region non- 
perturbative effects become very large, spoiling the evaluation of the perturbative part of the 
FF's: the PFF's evolved to NLO become indeed negative in the large x region. We have 
therefore preferred not to include that point in the fit. 

We have also presented along with the fits to ARGUS the results of a similar fit to OPAL 
data. Also in this case a few points have been excluded from the fit: the last one, where 
again large non-perturbative effects set it, and the first three ones, where the rise due to gluon 
splitting is observed. Since unaccounted for threshold effects may play an important role here, 
and the theoretical curve cannot be made to describe the data very well, we have preferred to 
avoid biasing the fitted parameters and therefore excluded this region altogether. 

The main result is the consistency of the two sets of parameters: the same values which 
fit the ARGUS data also describe the OPAL data, taken at a center-of-mass energy almost 
one order of magnitude larger. This finding lends support to our initial hypothesis of scale 
independence of the non-perturbative part of the fragmentation functions. 

Other fits with this "simple" non-perturbative FF have been performed, this time excluding 
the resummation of Sudakov terms. The reason for this is that when making convolutions of 



A. = 100 MeV 



A. = 200 MeV 



A. = 300 MeV 



Next-to-leading order evolution 



ARGUS 
OPAL 



.031 (1.09) 
.033 ± .005 (1.25) 



.019 (1.27) 
.015 ± .002 (1.54) 



.011 ± .003 
.008 ± .001 



1.53) 
1.72) 



Leading order evolution 
ARGUS .07 (1.65) .055 (2.1) 

OPAL .10 ± .01 (2.02) .08 ± .01 (2.48) 



.036 (2.72) 
.06 ± .01 (2.98) 



Table 2: Results for the fitting of the e parameter of the Peterson FF to AR- 
GUS and OPAL data, for three different values of A5 and with next-to leading 
order coefficient functions and NLO or LO evolution of the PFF's. Sudakov 
resummation is not included explicitly, and is therefore effectly reabsorbed 
into the non-perturbative FF. The number between the round brackets is the 
X^ per degree of freedom of each fit. 



the PFF's with the Sudakov included in the x space (rather than in Mellin moments space as 
we do now) the integration convergence is much more difficult. We have therefore chosen to 
incorporate the effect of the Sudakov resummation into the non-perturbative part, with the 
results given in the lower part of Table |1]. Once more, full consistency is found between the fits 
to ARGUS and to OPAL data. The results of these two fits are shown in figure [l|. 



3.2 Fits with the Peterson form 

Fits to the same ARGUS and OPAL data have also been performed using the Peterson form 
d^) as the non-perturbative part of the FF's. The fit is in this case a two- rather than a three- 
parameter one, namely the normalization and the e parameter only. Using NLO evolution and 
coefficient functions, but again no Sudakov resummation, and three different values for A5, we 
have found the results displayed in Table 0, while the curves resulting from these fits, for the 
choice A5 = 200 MeV, are shown in figure |I]. 

It is to be noted that the fitter was not able, in a few instances, to produce realistic errors 
when fitting ARGUS data, due to numerical inaccuracies resulting from the inverse Mellin 
transform of the Peterson FF. However, taking the error in the corresponding ffi to OPAL data 
as an indication, we see that also in this case the two fits are consistent, pointing to a scale 
independence of the non-perturbative part of the fragmentation functions. 

The most striking feature of these fits is however the discrepancy between their results and 
the value commonly used for the parameter e when describing c quarks fragmentation to D* 
mesons. It is indeed found in the literature (see, for instance, Ref. []2T|), and has been used 
in recent phenomenological papers [p2| , |T^, the value e = 0.06. The fitted values (except for 
the one at A5 = 100 MeV) also appear to be at variance with the result found by the OPAL 
Collaboration pO| as a fit to their own data, eopAL = 0.035 ± 0.007 ± 0.006. 



This discrepancy should however not come as a surprise if one considers carefully how e 
so far has been extracted from experimental data. Experiments usually report the energy or 
momentum fraction [xe or Xp) of the observed hadron with respect to the beam energy. On the 
other hand the momentum fraction which appears as the argument of the non-perturbative FF 
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NLO QCD + Peterson, e = 0.015 
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Aj = 200 MeV 
Normalization: fitted 




Figure 1: Distributions of D* mesons as measured by the ARGUS and OPAL 
experiments, together with the theoretical curves fitted to the same data with 



the (1 



a (3 



X) X' 



(full line) and the Peterson (dashed line) non-perturbative 



fragmentation functions. 



is rather the fraction with respect to the fragmenting quark momentum, usually denoted by z 



(see for instance ||21| for a discussion on this point). These two fractions are not coincident, due 
to radiation processes which lower the energy of the quark before it fragments into the hadron. 



In order to deconvolute these effects one usually runs a Monte Carlo simulation of the collision 
process at hand, including both the parton showers and the subsequent hadronization of the 
partons into the observable hadrons. The latter can be parametrized in the Monte Carlo by 
the same Peterson fragmentation function we have been using, and the value of e which best 
describes the data can be extracted. But what can be different in our approach is of course the 
perturbative QCD part, namely the parton shower. This showering softens the distribution of 
the partons, producing an effect qualitatively similar to that of the non-perturbative FF. On 
the quantitative level, the amount of softening (and hence the value of e) required by the non- 
perturbative FF to describe the data is related to the amount of softening already performed at 
the perturbative level. Monte Carlo's simulations so far only perform a leading order description 
of the showering, and can hence differ from our NLO evolution. 

Therefore there is not a "unique" and "true" value for the parameter e, but only a value 
closely interconnected with the details of the description of the pQCD part of the problem. 
For instance, a higher value of A5 results in a larger as and hence in more parton showering. 
This softens even more the perturbative part of the FF, and consequently less softening will be 
required from the non-perturbative part. The results in Table |^ show that this is indeed the 
case, a smaller value of e corresponding to a harder Peterson FF. 

A double check that the different description of the perturbative part can indeed responsible 
for the different e can be done by rerunning our fits with a leading order evolution, in such a 
way to mimick as closely as possible the Monte Carlo description of the process. The results 
are displayed in Table |^, and can be seen to be indeed much closer to the commonly used value 
of 0.06. The tendency to a discrepancy between ARGUS and OPAL fits could actually be an 
indication of the inadequacy of a leading order description of the scale violations taking place 
from 10 to 90 GeV. All this should however not be taken literally, as many other details might 
be included in the leading order Monte Carlo description of the perturbative showering and be 
missing or differently treated here. 

A further check of the modification of the value for e when going from a leading to a next- 
to-leading description of the perturbative parton shower can be obtained in the following way, 
to be taken as a kind of toy-model. 

Consider a distribution for the energy variable x, like the ones given by ARGUS and OPAL 
and plotted in fig. |I]. Thinking of them as described by the convolution of a perturbative and 
a non perturbative fragmentation function, the average value of x, call it {x)exp, can be written 
as a product of the average values of the perturbative and the non-perturbative FF's, i.e. 

If we now assume that both a leading and a next-to-leading description of the perturbative 
part can describe the data, provided they are matched by the appropriate non-perturbative FF 
(i.e., the appropriate value of e is chosen), we can write 



\X)exp \-^)pert\-^)np \-^ I pert \-^lnp i \^^) 



which leads us to 



K"^ I pert 



(r'^NLO _ \-^lpert , ,LO /. p-n 

\-^lnp l^\NLO\ Inp- X^-^ ) 



In this equation {x)pert refers to the second MeUin moment of the perturbative fragmentation 
function D^, while the {x)np can be calculated from the Peterson FF, like {x)np = J xD{x; e)dx. 
The suffixes "LO" and "NLO" on the perturbative parts mean that a leading or next-to-leading 
evolution kernel has been included before taking the average. The non-perturbative part is 
considered to be adjusted to fit the data together with the given perturbative term. 

The perturbative fragmentation function returns the following averages when evolved with 
A. = 200 MeV: 





\^ I pert 


Y^lpert 


\-^lpertl \'^l pert 


10.6 GeV 


.75 


.65 


1.15 


91.2 GeV 


.64 


.56 


1.14 



We can clearly see from this table how the NLO description does indeed soften the perturbative 
FF more than the LO one, producing a lower value for the average energy. 

Assuming e = 0.06 to be the right value to describe the data when a leading order pertur- 
bative description is used, we get {x)^^ = 0.67 and hence, from eq. (]T5|), {x)^^^ = 0.77. Upon 
inspection we see this average value for the Peterson FF corresponds to e = 0.016, i.e. a value 
fully compatible with the ones returned by the fits. 

Before closing this Section on the fits, we wish to point out once more that there is not a 
"best candidate" value for e, but only a value of e more suited to match the description of the 
perturbative showering one is actually employing. Surely enough, if the QCD description is at 
NLO a harder e, like our e = 0.015, should be used rather than the larger (and softer) e = 0.06, 
since part of the softening is now already included through more perturbative gluon emission. 



4 Production in ep collisions 

The use of fragmentation functions for heavy quarks to evaluate NLO cross sections for charm 
photoproduction has already been considered in Ref. 0. 

In this paper we use exactly the same formalism to evaluate cross sections for D* production, 
by complementing the PFF's used in the previous work with a non-perturbative component as 
described by eq. (|l|) and according to Ref. 0. 

The 7J9 cross section reads, schematically, 

da^p = / Fi/pda^i^kDj. + / Fi/pFj/^dcXij^kD^. . (16) 

In this expression Fi/p and Fj/^ are the parton distribution functions (pdf's) for the proton 
and the photon, since the so called direct and resolved component are both included. Unless 



otherwise stated, we will make use of the MRS-G |2^ and ACFGP |24| sets respectively. The a's 
are the kernel cross sections (= coefficient functions) for massless parton production ^, ^j and 
D^ is the meson fragmentation function of eq. (|l|). We will use in this FF the non-perturbative 
parameters fitted in the previous section to e~^e~ data and, since the non-perturbative FF's 
are normalized to one, we include the branching ratio BR{c — * D*^) = BR{c -^ D*^) = 0.26 



20|. This produces an absolute, parameter free, prediction, to be directly compared with the 



experimental data. 



We also convolute our 7p cross sections with the Weizsacker- Williams flux factor, 



0", 



ep\ 



Vmax 



s) = / dyf^/e{y)a^p{ys) (17) 



ymin 



with 

fi/e{y) 



a 

2^ 



^ + i^-y)\^QL. , o^2.. / 1 



In ^^ + 2miy 



y ^min \ ^max ^r. 



:i8) 



where y = E^/Ee, Qmin = "^eiz/ll ~ v) ^-nd mg is the electron mass, to mimick as closely as 
possible the experimental setup. For comparisons with ZEUS data we will adopt Q^ax = 4 
GeV^ and ymin = 0.147, ymax = 0.869, according to |Q. Moreover, we will present cross 



sections in the pseudorapidity range — 1.5 < r/ < 1 and in the pr range 3 < pt < 12 GeV. 

As already stressed in Ref. 0, it is important to point out how this low px boundary 
casts doubts on the validity of an approach based on the use of massless cross section kernels, 
and which had originally been devised for the resummation of large logarithms in the large pt 
region. In principle, one is missing terms of order m/pT, and the errors may therefore be large 
when Pt — rn. Only a comparison with a full massive calculation can finally assess whether 
the results are meaningful enough. Such a comparison will be presented in fig. §. 

A description of D* photoproduction in ep collisions similar to ours has recently been given 
Jl^]. When including the Peterson FF these authors tackle the problem from an apparently 



m 



different point of view, by evolving directly this non-perturbative FF and inserting instead the 
initial conditions @, (^ and (|^) for the heavy quark PFF's into the coefficient functions for 
7P to massless parton scattering. One can however easily see the two approaches are equivalent 
at the perturbative level. The Appendix does indeed show how they should only differ by 
uncontrollable higher order terms and, other than this, in the interpretation of the various 
components. 

Therefore the approach introduced in Ref. and now discussed here in detail, and the 
one successively used in |T^ should give similar results. We compare them in fig. |^. It shows 



the curve extracted from Ref. [T^ (wide-dotted line) and our results, for the same value of 
e = 0.06. No agreement is found, however, neither (full line) with what will be our standard 



choice of renormalization/factorization scales {fi = fi^ = fip = ^"^t = ^\lm? + p^ with ^ = 1 
/io = TTt, and A5 = 200 MeV) nor (dashed line) when we make the same choice as Ref. [p!7[ 
taking /^/j = mT-, /i^ = Sm^, /xq = 2m, GRV-G HO ||2^ as the photon pdf's set. 



Spurious higher order terms could be responsible for the discrepancy. If one does indeed 
check fig. 2 of Ref. |T^, by comparing curves C and D a difference similar to the one found 



above between the wide-dotted and the dashed line can be seen. This large difference could 
therefore be due to the moving of the initial condition terms for the fragmentation function to 
the kernel cross sections for massless parton scattering (see Appendix). Curve D of Ref. [117 



has been made following our standard PFF formalism, and by comparing it with our results 
we have indeed found agreement. 

It is worth noting that the spurious terms contain large Sudakov logarithms of the form 
log(l — x), and could indeed be not negligible. Since we fitted e"'"e~ data with the same overall 
fragmentation function we are now using here, we believe the large effect of these terms - if 
present - to be effectively absorbed into the fixed non-perturbative component. Hence it should 
not spoil a reliable evaluation of photoproduction cross sections. 

10 



^9 

D 











ep, (27.5 + 820) GeV 


standard, ^ = 1 


- 


3GeV<p^<12GeV 


GRV-GHO,^g=l, 




ZEUS untagged setup 


^P = 2,Ho = 2m 






from KKS paper 


- 


MRS-G and ACFGP 


full massive calc. 


- 


e = 0.06 




-. 






\ 


^"^ 


~^^^^ 









0.75 1 



Figure 2: Comparison of our results with those of Ref. |T^ (KKS) and 



with the full NLO massive calculation of Ref. |2|]. The GRV-G HO photon 
parton distribution functions set is employed for all the curves except for the 
"standard" one (full line), ^/j and ^p refer to the ratios of the renormalization 
and factorization scales to the transverse mass m^ respectively. 



Also shown on the same plot is the result of the full NLO massive calculation by |^ (close- 
dotted line), itself convoluted with a Peterson FF with e = 0.06 too, as taken from ||2^. Good 
agreement with our result is found, especially when making our standard choice of scales. Such 
a successful comparison could probably not have been expected beforehand, given the missing 
m/pT terms, but a posteriori it can perhaps be considered a check of our results, being the 
massive result the benchmark at these low pt values. The agreement will also allow us to 
extrapolate to the massive calculation the effect of varying the value of e. 



4.1 Comparison with experiment 

We now compare our results with experimental photoproduction data obtained at HERA by 
ZEUS |26[ and HI ^ Collaborations. 

We first plot, in fig. |], the pseudorapidity distributions, integrated over the pxi obtained 
with the Peterson FF with different values of e. These results have been obtained with the 
pdf's set MRS-G for the proton, and ACFGP for the photon in the resolved component. For 
the renormalization/factorization scales we made the standard choice /i = fiR = fip = ttit and 
taken /iq = Tn as the starting value for the evolution of the FF's. A5 is taken equal to 200 MeV. 

As expected, the use of a smaller e hardens the non-perturbative FF and hence enhances 
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Figure 3: Pseudorapidity distribution of D* as measured by the ZEUS ex- 
periment and theoretical predictions for different values of e in the Peterson 
fragmentation function (full line) and with the (1 — x)°'x^ FF (dashed line). 



the cross section, since the partonic kernels fall rapidly with increasing pt- The cross section 
obtained with e = 0.015 is 50% larger than that with e = 0.06, and while the latter seems to 
fall short of describing the ZEUS data, the former does a good job, at least in the first two 
bins. But we emphasize here once more how a full assessment of the reliability of these results 
needs a comparison with the full massive calculation, rather than with the experimental data, 
which however need to be improved in precision. 

For comparison, the cross section obtained with the simple FF, (1 — x)°'x^, with a = 0.9 
and (3 = 6.4, is also shown (dashed line) in fig. |^. These values for a and (3 fit the OPAL data 
from e~^e~ collisions like e = 0.015 does, see Section |^, and the photoproduction cross sections 
are indeed also in good agreement. This on one side shows how in this case there is little 
dependency on the precise shape of the non-perturbative fragmentation function. On the other 
side, it strengthens our trust of the cross section with the Peterson, much harder to evaluate 
due to the numerical difficulties related to the inverse Mellin transform. 



The total cross sections, obtained by integrating the curves in fig. |^ over the pseudorapidity, 
are also shown in Table ^ They are to be compared with the experimental result from ZEUS 
26[ o" = 10.6±1.7{stat.)±\-^{syst.) nb. Notice that the 17% increase found going from e = 0.06 
to 0.035 is in good agreement with the 15% estimated in ||2H] using the massive calculation. 



To get a feeling of the stability of our results we plot in fig. |a the results obtained for the 
pseudorapidity distribution with different choices of renormalization/factorization scales and 
with the conservative value e = 0.02. While the central curve is obtained with /i = rriT, the 



12 





a (nb) 


a = 0.9, /? = 6.4 


8.0 


e = 0.015 


7.7 


e = 0.02, ^ = 0.5 


9.6 


e = 0.02 


7.2 


e = 0.02, i = 2 


6.2 


e = 0.035 


6.2 


e = 0.06 


5.3 



Table 3: Predictions for the total cross sections in the ZEUS untagged setup, 
3 < J9T < 12 GeV and -1.5 < ?? < 1. 



two others are produced with /i = C,mT, with ^ = 0.5 and 2. The variation is not negligible, 
especially in the lower scale direction, but we should bear in mind that at such a low scale we 
are at the border of the applicability of perturbative QCD. Also shown on this plot are the 
results obtained with two other photon pdf's sets, namely GRV-G HO and AFG [^]. The 
variations are smaller than those given by varying the scales. 

By comparing with the experimental results we can see that we can get a fairly good 
description of the data already with a central choice of scales. 

A similar comparison is also made, in fig. ^d, with the px distribution obtained by the 
ZEUS Collaboration. The curves, obtained with e = 0.02, seem to offer a fair description of 
the data. 

Finally, we want to present a comparison of the results of our approach with more sets of 
data. We now use HI results, both in the tagged and the untagged experimental setup. These 
we reproduce by taking in the Weizsacker- Williams convolution Ql^^x = 0-01 GeV^, ymin = 0.28, 
Umax = 0.65 and Q^^^ = 4 GeV^ ymin = 0.1, Vmax = 0.8 respectively, according to Ref. [^. 
Fig. H shows the results for the rapidity distribution^ , obtained with the Peterson FF with 
e = 0.02 and the standard choice of scales, i.e. ^ = 1. 

The total cross sections for these curves, integrated within the 2.5 < pt < 10 GeV and 
— 1.5 < y < 1 range, read 4.2 nb and 14.4 nb for the tagged and the untagged sample respec- 
tively, to be compared with the experimental results 4.9 ± 0.7±o;5g nb and 20.2 ± 3.3±3;g nb. 
A quite good agreement can be seen, especially for the tagged sample. 



5 Conclusions 



In this paper we have applied the technique of fragmentation functions for heavy mesons to D* 
production in e~^e~ and ep collisions. 

These fragmentation functions are made of a perturbative part, which we evolve with next- 



^ HI presents its experimental results as a function of the rapidity rather than of the pseudorapidity. Our 
approach, in that it deals with massless partons in the kernel cross sections, cannot distinguish between the 



two. The two quantities become of course identical in the large px region, and at pt 
differ by about 10%. 



2.5 GeV already only 
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malization and factorization scales, as fi = ^rriT, and of the photon pdf's 
sets. 
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Figure 5: Comparison of our results with the experimental data from the HI 
Collaboration, Ref. |2^]. 



to-leading accuracy, and a non-perturbative one, which we fit to e^e~ data taken by ARGUS 
and OPAL and subsequently use to predict photopro duct ion cross sections, to be compared 
with data by HI and ZEUS. 
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When fitting e+e^ data with a Peterson non-perturbative form we find values for the e 
parameter sensibly different from the commonly accepted value 0.06. A central value for our 
fits, when using NLO evolution, is e = 0.02. This hardens the non-perturbative fragmentation 
function, and increases the photoproduction cross section, bringing it in better agreement with 
the data. 

Our photoproduction results are found in good agreement with the NLO full massive ones, 
which are reliable at the low values of px probed by the experiments and can be taken as a 
benchmark for comparisons. Convoluting them with a Peterson with a lower e will also increase 
the cross section, again producing a better agreement with the data. Slightly less conservative 
choices than those made here for the renormalization/factorization scales, the photon parton 
distribution functions set, the c ^ D* branching ratio and the value of e could easily make the 
agreement even better. 
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A Appendix 

In this Appendix we show how the approaches of Refs. [Q] and |T^ are identical at the pertur- 
bative level. 

Consider a cross section for producing a heavy quark of mass m at the large scale Q, (j{Q, m), 
given by the convolution of a coefficient function C{Q,fi) and a perturbative fragmentation 
function D{fi,m). In the Mellin moments space we write this as a product: 

a{Q,m) = C{Q,fi)D{fi,m), (19) 

and /i is the factorization scale. Since D{fi, m) is the fragmentation function evolved up to the 
scale /i, we can write it in terms of an initial condition at a scale /iq as 

D{n, m) = E{jj, jjo)D{no, m)Dnp. (20) 

The factor £'(/i, /xq) is the so called evolution kernel, and we have now also included a non- 
perturbative term Dnp, for instance the Peterson FF, according to eq. (|I|). Indeed, to think it 
to multiply the perturbative initial condition or the evolved PFF is absolutely identical, since 
D(/i) is in both cases simply a product of three terms. 

Putting together the two equations we have 

cr(Q, m) = C{Q, /i)-E(/i, /io)-D(/io, m)Dnp, (21) 

which is for instance the way we write our e~^e~ cross section in Mellin space, the one in x-space 
to be found by numerical inverse Mellin transform. 

If we now consider that both the coefficient functions (see for instance Ref. |T^) and the 
initial conditions of the PFF's (see eqs. (§), (0) and (|^)) can be calculated as series expansions 
in as, like 

C((5,/i) = 1 + as(/i)c((5,/i) and D{fio,m) = 1 + as{iJo)d{fio,m), (22) 

inserting these expressions into eq.(^) and rearranging it, up to uncontrollable 0{a'^) terms 
we can write 

a{Q, m) = (l + as{fi)c{Q, /i) + as(yUo)rf(/io, 'rn)^E{j2, /io)-Dnp- (23) 

This is (with the exception of as{fio) which they take as{fi) instead) the form employed in 
Ref. |]T^ when E{fi, fio)Dnp is considered as an "evolved" non-perturbative FF, and with the 
(i(/io,'m.) functions changing the coefficient function's scheme. If one takes /iq = /i the new 
coefficient function will be close to the cross section for massive quark production, containing 
the logarithmic terms log(Q/m). Indeed, the d{^Q,m) functions had been determined in [Q 
exactly this way, but going the opposite way round, i.e. evaluating the full massive a{Q,m), 
extracting from it the coefficient function c{Q, /i) in the MS scheme, and defining the remaining 
piece to be the initial state condition of the fragmentation function. 
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